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We present a simple model of Majorana fermions on a square lattice, and study zero- 
energy states due to Z2 vortices. We show the relationship between the Chern number of 
the ground state and the number of zero-energy states through numerical calculations 
for finite systems. We also discuss the same relationship for the bulk system by observing 
the change in the spectral asymmetry as well as the lattice index. 
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Quantum statistics of identical particles is one of the most basic properties in the 
study of many-body systems. Among some unusual statistics of elementary excitations 
known so far, Majorana fermions obeying non-Abelian statistics are of much current 
interest and are expected to play a crucial role in fault-tolerant quantum computa- 
tions. ^"^^ Since the field theory of Majorana fermions^^ describes the critical point of 
the Ising model, they are often called Ising fermions. Kitaev proposed a solvable spin 
model on a honeycomb lattice, in which spins can be mapped to Majorana fermions. 
A Px + ipy superconductor is another well-known example of them.^'^^ These models be- 
long to the same universality class D for topological superconductors,^' and conversely, 
superconductors in class D may be regarded as Majorana fermion systems.^' 

Defined on a lattice, the Kitaev model provides not only various numerical results^^^ 
but also intriguing notions such as the Majorana number, which relates the existence 
of the zero-energy states in the Z2 vortex sectors to the Chern number of the ground 
state. 2) The purpose of this paper is to present a simpler minimal lattice model written 
directly by Majorana fermions and defined on a conventional square lattice. This model 
is intimately related to the Wilson-Dirac operator in the lattice gauge theories. The 
development in this field concerning a chiral anomaly on a lattice^^' yields a notion 
of the spectral flow and resultant Z2 index theorem for the zero-energy state of the 
Wilson-Dirac Hamiltonian. 
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Hamiltonian density is defined by 

-H = + v;) + 73 + ^v;v,^ , (1) 

where // = 1, 2 and the lattice derivatives are defined by 

aV ij,il){x) = Ui^{x)iIj{x + afi) - ipix), 

aV^^il^ix) ^ ijj(x) - U*(x - all)'ijj(x - ajl) (2) 

with the lattice constant a and the unit vector jl in the n direction. The link variable 

U^{x) is defined on the link between x and x + afi. Let us choose the 7-matrices as 
7i = o-"^, 72 = 0"^, and 73 = —a^, and restrict the link variables as Un{x) = ±1. Then, it 
turns out that 

n* = -n, (3) 

which means that the Hamiltonian is purely imaginary and antisymmetric. This enables 
us to define the following Majorana fermion model on the square lattice: 

X 

j 

where ip{x) in the first line denotes a Majorana fermion operator on the site x obeying 
{ip{x),ip{y)} = Sxy/a. In the next lines, we have introduced the dimensionless operator 
Cj = a^^'^^(x) as well as dimensionless parameters M = ma and B = b/a, where j = x/a 
is the set of integers indicating the lattice point. 

It may be noted that the Hamiltonian density (1) is related to, as will be 
discussed momentarily, the Wilson-Dirac operator on the square lattice, Dw — 
[7^4 (V^ + V*) — 6V* V^] /2, where the Laplacian term is the famous Wilson term in- 
troduced for the purpose of avoiding species doubling. Remarkably, the Hamiltonian (4) 
is intimately related to condensed matter physics, which has been used as a prototype 
of topological insulators. ^^"^^^ Although all these studies focus on the Dirac fermions 
on a lattice, the Hamiltonian density (1) serves as a simple Majorana fermion model, 
similarly to the Kitaev model. 

Let us begin with discussions on the Chern number of the ground state without 
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vortices. To this end, we set U^{x) — 1 for all x. The Fourier transformation yields 



73 = Xa'Ja, (5) 



M L M - 

where and {fi = 1, 2) are lattice momenta defined by = sin /c^ and = cos /c^, 
and A runs A = 1, 2, 3. The energy spectrum is given by ±R, which follows from Hi — 
with R — a/^^ X\. The first Chern number c characterizing the negative energy 
state of eq. (5) may be computed using the projection operator = (1 — Hk/R)/'^ — 
(1 — Xa'Ja) /"^ such that 

c=±ltrPiiPr^lJn, (6) 

where Xa is defined by Xa = Xa/R, which is regarded as the coordinates of because 
of — l.The exterior derivative d is with respect to the momentum defined by 

d = dk/j^dk^ , and Q is a two- form given by 

^ = -^^ABcXAdXBdXc 



C1C2 



d^k. (7) 



As discussed by Fujiwara et a/./^^ k^ Xa defines the mapping / : T^ — )■ S^, and 
the integral of Vt over the Brillouin zone in eq. (6) is simply the degree of the mapping 
/ times the volume of -S"^, i.e., 47r. To calculate the degree of the mapping /, let us 
consider the point X — (0, 0, 1) on S^. The preimage of this point under / is the set of 
points k^ = 0, TT. Note that X^ = sgn(M) for k = (0,0), sgn(M - 2B) for (0,7r) and 
(tt, 0), and sgn(M — AB) for (tt, tt). Taking the sign of the prefactor C1C2 in eq. (7) into 
account, we finally conclude that the Chern number is c = sgn(S) forO < M/{2B) < 1, 
c = -sgn(i?) for 1 < M/{2B) < 2, and c = otherwise. 

We now begin the discussions on the zero-energy states around Z2 vortices. For 
numerical calculations of the finite size systems, we set two Z2 vortices separated by 
distance d, as illustrated in Fig. 1(a). Since the spectrum is symmetric with respect 
to the zero-energy owing to property (3), we denote the eigenenergies as ±e„ {n — 
0, 1, 2, • • • ). Figure 2(a) shows the spectrum of the system that has a,c—l ground state. 
There appear two approximately zero-energy states in the presence of well-separated 
Z2 vortices. However, for smaller d, these states have finite energies, ±eo, owing to their 
interactions. We show in Fig. 2(b) the energy ±eo as a function of the distance d. The 



J. Phys. Soc. Jpn. 



LETTERS 



- -X 



d 



(a) 



(b) 



Fig. 1. (a) Two Z2 vortices separated by distance d. On the thick links crossing the dashed line, 



Uu. 



-1 link variables are set. The crosses represent the positions of Z2 vortexes, (b) Z2 link variable 



with = —1 in (a) is replaced by a generic U(l) variable e*"^. Because of the complex link variables, 
an arrow indicating the positive (j) direction is attached to each link, and the curved arrows show the 
positions of U(l) vortexes. 



energy cq can be well parameterized as^^'' 

eo = Ae-^/«, (8) 

where A is of the order of the bulk energy gap and ^ is the coherence length between 
two Z2 vortices. This behavior is in sharp contrast to the Kitaev model, which shows 
an oscillation of A as a function of d}^^ Since the coherence length generically depends 
on the model parameters, we plot ^ as a function of M in Fig. 2(c). 

We next show the zero-energy states from the point of view of spectral flow, by 
extending the Z2 link variable with = —1 to the generic U(l) variable = e*"^. 
As illustrated in Fig. 1(b), this yields a pair of a U(l) vortex and an antivortex, and 
exactly at = vr, they become the two Z2 vortices discussed above. We show in Fig. 
2(d) the spectral flow as a function of (p. Starting at = 0, which corresponds to the 
ground state without vortices, and increasing adiabatically, two eigenvalues approach 
zero, and at = vr, they form exponentially degenerate zero-energy states. The in- 
creasing (decreasing) state denoted by red (blue) points is associated with the vortex 
(antivortex). This behavior suggests that if there is a single vortex in the infinite bulk 
system, only one state (in the present case, only the increasing state denoted by red 
points) crosses the zero energy without a level repulsion, and hence, we have just one 
zero-energy state sitting exactly on the zero energy at = vr. 

We have so far shown some results of the model with a c = 1 ground state. We 
mention in passing that contrary to this case, we numerically confirm that the model 
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Fig. 2. (Color online) (a) Spectrum of finite size system (20x 20) under periodic boundary conditions. 
The distance between vortices is d = 9, and other parameters used are M = 1 and _B = 1 (c = f for 
ground state). Inset: Spectrum around the zero energy, (b) Energy ±eo as a function of d. (c) Coherence 
length ^ as a function of M. (d) Spectral flow as a function of <j) when all Z2 link variables with t/^ = — 1 
are replaced by Uf^ — e*"^. 



with a c = ground state shows no zero-energy states even in the presence of distant 
vortices. 

We now elucidate a true zero-energy state associated with a single Z2 vortex in the 
bulk system. We here remind ourselves that the Hamiltonian density eq. (1) is intimately 
related to the Wilson-Dirac operator. From the point of view of lattice fermions, this 
operator is free from the species doubling owing to the Wilson term, but at the cost of 
a loss of chiral symmetry. Recent developments, however, enable us to define a variant 
of chiral symmetry on a lattice. A key property of such developments is the Ginsparg- 
Wilson (GW) relation, which has been recently realized as the overlap Dirac operator. 

The GW relation is given, in the case of two dimensions, by 73D + D73 = aD'-f^D, 
where D is a generic Dirac operator obtained by the renormalization group for a lattice 
fermion.^^^ In this equation, the r.h.s means broken chiral symmetry since the lattice 
constant is finite. Suppose that D satisfies the GW relation and that 73 is a variant of 
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73 on a lattice defined by 73 = 73(1 — aD). Then, one can show that 73D + D% — 
owing to the GW relation. This can be regarded as chiral symmetry on a lattice. A 
concrete operator satisfying the GW relation was constructed and called the overlap 
Dirac operator. It is defined by the use of the Wilson-Dirac operator as follows. Let 
A be an operator defined hj A = M — aDw- Then, the overlap operator D is defined 
hy D = [1 — ^(^"1"^)"^/^] /a. With the use of the chiral anomaly can be reproduced 
on a lattice such that q{x) — tr [73 (1 — aD{x, x)/2)] /a^, which also leads to the index 
theorem on a lattice.^^^ 

With appropriate choices of the 7-matrices for T-L in eq. (1) and as well, we have 
oH — 73A, and hence the overlap Dirac operator can be expressed as 

D = ^(l- 'ysH/Vw) . (9) 

The chiral anomaly can then be expressed by q{x) — tr ^73 + T-L/VT-{P^ /(2a^)- This 
has a topological meaning through the index theorem mentioned above. This motivates 
us to study 

Q^a'J2<li^)-lT^^^ (10) 

where "H is the Hamiltonian density defined by eq. (1), and in the last equality the 73- 
term independent of T-L has been omitted, since it vanishes. Alternatively, this quantity 
can be viewed as a spectral asymmetry, Q — Tr (P+—P_)/2, where P± = {l±H/Vy?)/2 
represents the projection operator to the positive or negative energy state. 

To explore the zero-energy state around a vortex, we set a single Z2 vortex at the 
origin accompanied by an infinite sequence of — 1 links from the origin to infinity. In this 
case, eq. (10) is ill-defined, since H is expected to have a zero-energy state. To overcome 
the difficulty, we utilize spectral fiow again. Namely, we extend the link variables with 
U^{x) = — 1 to a generic U(l) variable, U^{x) = e*"^, to see spectral fiow as a function of 
0. This idea was originated by Roy,^^-* who studied the zero modes on the basis of the 
arguments of Laughlin^^^ and Halperin.^^^ In what follows, we will give an alternative 
proof of spectral flow, by calculating it directly. 

The extended Hamiltonian density, denoted by 1-L{(p), becomes non-antisymmetric, 
and property (3) is modified as 

?^*((/.) = -n{-(i>). (11) 

For special cases = 0, tt, and 27r, the extended Hamiltonian describes the Majorana 
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fermions: 7^(0) and 'H{2n) are the original vortex-free Hamiltonian, and Hin) is the 
Hamiltonian in the presence of a Z2 vortex. Of course, T-L{0) and 'H{2n) are identical. 
However, if the spectrum is followed as a function of 0, spectral flow occurs and the 
states are generically rearranged. Even when "H has zero eigenvalues at some 0, in 
particular, at = tt, we can avoid them by adding some infinitesimal perturbations. 
For example, around (j) — n, we extend into — )■ + ie, and then, the Hamiltonian 
density becomes non-hermitian and spectral flow bypasses the zero eigenvalues. Thus, 
we can compute the spectral asymmetry (10) continuously from = Oto0 = 27r — 0. 

Of course, it is difficult to calculate it directly for Hamiltonian (1). Remarkably, its 
general form was determined by Liischer^^^ using only the cohomological arguments on 
the basis of the locality, the gauge invariance, and the topological nature of the anomaly. 
Some unknown parameters in this generic form were determined by direct calculation 
up to a°, assuming a smoothly interpolated gauge potential on the links. In the plane 
wave basis, q{x) of 0{aP) can be computed from 

q{x) = - I - — -rtr — , (12 

^ ' 2 7 (27r)2 ^/WT{HiJW+^ 

where the leading term T-Lk is defined by cq. (5), and dl-L is the remaining difference 

operators defined by the relation er'^^^^'He^^^ — Hk + SH. Note that tr 7^7^7(7 = 2ieABC: 

the leading contribution in eq. (12) is given by 

I r dPk 1 



^^F,, + 0{a% (13) 

where c is the Chern number in cq. (6), and F12 is the field strength associated with the 
U(l) gauge potential. It should be stressed here that the above F12 is the field strength 
of order a°, i.e., in the continuum limit, assuming that the U(l) link variable can be 
parameterized by a smooth gauge potential. Having a known leading order, we obtain 
the exact g(x), according to the Liischer's arguments, as 

q{x) = ^F,2{x) + d;k,{x), (14) 

where kn{x) is a gauge-invariant local current, 9* is V* with C/^ = 1 in eq. (2), and F-12 
is the field strength on a lattice: 

a^Fi2{x) = -ilnUi{x)U2ix + ai)U*{x + a2)U;{x). (15) 

It should be noted that the current k^ is a lattice artifact,^^) and vanishes in the con- 
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tinuum limit a — > 0. 

Now, let us explore the zero mode in the presence of a single Z2 vortex. To this 
end, we investigate the spectral flow as a function of 0. First, when = 0, the model 
includes no vortices. We readily sec F12 = 0. Note also the fact that the Hamiltonian has 
a well-defined continuum limit, implying that o^^^dfj^k*^ = 0. Therefore, we conclude 
that Q = when = 0. Let us now increase adiabatically from = 0. We specify the 
branch of the logarithm in eq. (15) such that < a^Fi2 < 27r. Then, the field strength 
(15) becomes a^Fi2 = at the plaquette of the vortex and a^Fi2 = elsewhere. 
Unfortunately, for a finite 0, it is impossible to take the continuum limit of the model 
for which it is difficult to calculate and therefore, Q. Finally, when — > 27r — 0, it 
follows from eq. (14) that Q ^ c, since in this limit, the continuum limit is well defined 
and the contribution from the current A;^ should vanish again. Therefore, we obtain the 
spectral flow 

5Q = g(27r - 0) - Q(0) = c. (16) 

More generically, if one considers the case of n± vortices with ± vorticity, SQ = cSn, 
where dn = — n_. 

For simplicity, we restrict our discussions to the case of c = 1. Then, 6Q = 1, 
which implies that one negative state at = flows towards the positive energy as 
increases, and flnally at = 27r, it becomes a positive energy state above the gap. Let 
e(0) be the eigenenergy of such a state. Then, the symmetry (11) guarantees that the 
spectral flow is symmetric with respect to = tt, e(0) = — e(27r — 0), and therefore, 
e(7r) = 0. This state is the zero-energy state due to a Z2 vortex. On the other hand, 
if the Chern number of the ground state is trivial, c — 0, the spectral flow does not 
occur and no zero-energy state appears. All these observations are consistent with the 
numerical calculations. More generically, if SQ = even, an even number of states flow 
and they cross the zero energy at generic owing to the level repulsion at = tt, 
whereas if 6Q = odd, generically, only one state should cross the zero energy exactly at 
= TT. Therefore, the number of zero-energy states associated with Z2 vortices is given 
hy SQ modulo 2. 

Although we have investigated eq. (10) as the spectral asymmetry, it should also 
have the meaning of the index of the Hamiltonian, and eq. (16) can be interpreted 
as the topological change of the index at = tt. To explore the index, we note that 
only the zero-energy states of the Dirac operator D in eq. (9) contribute to the index 
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(10).^^^ In what follows, for simplicity, we restrict our discussions to the nontrivial 
case < M < 2B, giving c 1 for the ground state. Then, only one mode around 
k ~ (0,0) becomes massless in cq. (9), whereas other doublers k ~ (0,7r), (7r,0), (7r,7r) 
are still massive. Therefore, the massless mode around zero momentum is expected to 
contribute to the index. To compute the index, we must introduce the flux for this 
mode. In the lattice model, the flux is located at one plaquette, which makes it difficult 
to consider the continuum limit of the model. Instead, we consider a flux defined in a 
wider region on a two dimensional space, which should give the same index as long as 
the total flux is kept unchanged. Thus, we explore the effective continuum Hamiltonian 
that we expect gives the index (10), 

^c = -^7/.(5/.-^^/.), (17) 
where we assume that the gauge potential is given by^^^ 

A^{x) = e^,x,A{r) (18) 

with the asymptotic form 

f (r^O) 
A{r)^{ (19) 

Integrated over 2D, this gauge potential indeed yields the total flux 4>. 

Suitable unitary transformation enables us to choose 7i = cr^, 72 = cr^, and 73 = o"^. 
In this representation, let us obtain the zero-energy states T-Lci^ — 0, which is explicitly 
given by each chiral component ± such that 

dr±-deTA{r)]iP±{r,e) = 0, (20) 
r / 

where ip±{r, 9) is the zero-energy wavefunction with chirality ±, i.e., ^zi^± — i^'i- Set 
'0±('^, — e^^"^^'ip±{r). Then, the equation is reduced to 

(dr-jTAir))^±ir) = Q. (21) 

We obtain the solutions of this equation with the following power law behavior at the 
boundaries, 

I (r ^ 0) 

V'± ~ < • 22 

I ^m=F0/(27r) _^ 

The normalizability of the wavefunction requires 

-1< m < ±(/'/(27r) - 1. (23) 
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When = 0, there are no solutions, but when exceeds tt, there appears the solu- 
tion m — —1/2 in the chirality + sector. This wavefunction is antiperiodic (double- 
valued),^'''^ ijj+iO = 2%) = —il^+iO = 0). The change in its index from to 1 (16) owing 
to spectral flow may thus be interpreted by this effective model. Although the model 
(17) corresponds to the Dirac fermion with U(l) gauge potential, it can be written as 
Majorana doublet fermions with an 0(2) gauge potential. 
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